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l—j , I calculate the diffuse flux of electron antineutrinos from all supernovae using the 

' information on the neutrino spectrum from SN1987A and the information on the 

rate of supernovae from direct supernova observations. The interval of flux allowed 
at 99% confidence level is ~ 0.05 — 0.35 cm~ 2 s~ 1 above the SuperKamiokande (SK) 
^f) ' energy cut of 19.3 MeV. This result is at least a factor of ~ 4 smaller than the 

current SK upper limit of 1.2 cm _2 s _1 , thus motivating the experimental efforts to 
, lower the detection energy threshold or to upgrade to higher volumes. A Megaton 

water Cherenkov detector with ~ 90% efficiency would record ~ 2 — 44 inverse beta 
^ ■ decay events a year depending on the energy cut. 
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1 Introduction 

> ■ 

X: 

Neutrinos from core collapse supernovae are unique messengers of information 
on the physics of supernovae and on the properties of neutrinos. About the 
former, neutrinos are precious to study events that occur near the core of the 
star, where matter is opaque to photons: the neutronization due to electron 
capture, the infall phase, the formation and propagation of the shockwave and 
the cooling phase. Moreover, they allow to test the cosmological rate of super- 
nova neutrino bursts and thus to probe indirectly the history of star formation. 
Within neutrino physics, one can learn about the hierarchy (ordering) of the 
neutrino mass spectrum, about the e-3 entry of the neutrino mixing matrix, 
about possible non-standard neutrino interactions, existence of sterile states, 
etc.. 

The experimental study of supernova neutrinos is challenging in many re- 
spects. With current and upcoming neutrino telescopes, two scenarios are 
possible. The first is the detection of a burst from an individual supernova 
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which is close enough to us to produce a significant number of events in a 
detector. This limits the candidate stars to those within few hundreds of kilo- 
parsecs from the Earth, where the rate of core collapse is as low as ~ 2 — 3 
per century (see e.g. [1,2]). This explains why only one neutrino signal of 
this type has been recorded so far, in 1987 from SN1987A [3,4,5]. The second 
possibility is to study the diffuse flux of neutrinos from all supernovae. This 
requires massive detectors with a difficult background rejection [6,7] or very 
precise geochemical tests [8,9]. So far, the searches for this flux have given 
negative results, and upper limits were put. Among them, the most stringent 
is given by SuperKamiokande (SK) [6] on the flux of electron antineutrinos 
above E = 19.3 MeV in neutrino energy (see also the looser constraint from 
KamLAND in the energy interval 8.3 - 14.8 MeV [7]). This limit, at 90% 
confidence level, is: 

$(£ > 19.3 MeV) < 1.2 cnrV 1 , (1) 

and holds for a variety of theoretical inputs [6] 1 . The bound (1) approaches 
the range of theoretical predictions, and thus motivates the expectation that 
a positive signal may be seen in the near future, either with more statistics 
at SK or at the next generation Cherenkov detectors with Megaton volumes 
(20 times larger than SK) like UNO [10,11], HyperKamiokande (HK) [12], and 
MEMPHYS [13]. 

When assessing the possibility of a future measurement, one should consider 
that predictions of the diffuse supernova neutrino flux (DSNi/F) suffer large 
uncertainties, due to our poor knowledge of the underlying physics. In partic- 
ular, there are two sources of theoretical error. One is the uncertain value of 
the supernova rate 2 (SNR for brevity) as a function of the redshift. This rate 
can be obtained directly from supernova observations [14,15,16], or inferred 
from the star formation rate (SFR), which in turn can be extracted from the 
data on optical or far ultraviolet luminosity of galaxies [17,18,19,20]. It can 
also be constrained from the metal abundances in our local universe (see e.g. 
the discussion in [21]). These different methods have uncertainties of various 
nature: statistical, systematic, or due to theoretical priors. In particular, the 
connection between the SFR and the SNR relies on a number of theoretical 
inputs, such as the minimum mass required for a star to become a supernova. 
These inputs are uncertain and thus contribute to the error on the SNR, and 
ultimately to the error on the DSNz/F. 

1 The bound in Eq. (1) depends on the energy spectrum of the neutrinos arriving at 
Earth. Thus it is model-dependent, even though different models of spectra happen 
to give similar values of it [6] . 

2 In this context supernova rate means the rate of neutrino-emitting supernovae. It 
does not include objects that do not produce neutrinos in significant amount, such 
as the type la supernovae or core collapse supernovae that evolve into black holes 
before emitting neutrinos. 
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The second source of error on the DSNz/F is the uncertainty on the neutrino 
fluxes originally produced inside a supernova. To reduce this, one can rely 
on the results of numerical calculations of neutrino transport. An alternative 
possibility is to use the only experimental information available, i.e. the spectra 
of the neutrino events from SN1987A, as was first proposed by Fukugita and 
Kawasaki [22]. 

Motivated by theory or by observational results, several authors have combined 
representative neutrino spectra with realistic models of the SNR to obtain the 
DSNz/F [23,24,25,26]. Others have investigated the connection with neutrino 
detection [27,28,29,30,31,32,33] or the possibility to constrain the SNR and the 
SFR using the bound (1) [22,21]. In most of the calculations of the DSNz/F 
available in literature the SNR is inferred from the SFR and the neutrino 
spectra from numerical calculations are adopted. Exceptions are ref. [22], and 
ref. [26], where the SNR was constrained using the metal enrichment history 
of the universe. In all previous works, the quoted uncertainties on the DSNz/F 
are indicative. 



In this work I develop the complementary method of Fukugita and Kawasaki. 
Specifically, here I analyze the SN1987A data, taking into account neutrino 
oscillations, to constrain the neutrino fluxes in the different flavors produced 
inside the star. I also perform a fit of the measurements of the SNR from direct 
supernova observations. This is meant to be a first step towards a combined 
fit of all the available data, and has a value of its own because it is free from 
the uncertainties that effect the SFR-SNR connection, as emphasized in [21]. 
Finally, I combine the results of the two analyses, and use them to find the 
interval of values of the DSNz/F allowed at a given confidence level. This study 
answers the well defined question of what we can conclude on the DSNz/F if 
we decide to rely solely on direct experimental information. If compared to 
the previous literature, it shows how the prediction of the DSNz/F changes 
with a change of approach in the calculation, and this is important to provide 
robust guidance for experimental searches of this flux. At a more technical 
level, this work provides a statistically meaningful error on the DSNz/F. This 
is relevant considering that, as the technology of neutrino telescopes advances, 
likely the phase of discovery of a supernova neutrino signal will be replaced 
by a phase of detailed analyses and the consideration of uncertainties on the 
DSNz/F predictions will become necessary. 

The text is organized as follows: after a section on generalities (Sec. 2), I 
present the analysis of the SN1987A data in Sec. 3, and that of the SNR 
measurements in Sec. 4. The combination of the two and the final results for 
the DSNz/F are given in Sec. 5. Discussion and conclusions follow in Sec. 6. 
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2 Generalities 

Core collapse supernovae are the only site in the universe today where the 
matter density is large enough to have the buildup of a thermal gas of neutri- 
nos. Thanks to their lack of electromagnetic interaction, these neutrinos can 
diffuse out of the star over a time scale of few seconds, much shorter than the 
diffusion time of photons. This makes the neutrinos the principal channel of 
emission of the O(10 53 ) ergs of gravitational energy that is liberated in the 
collapse. The energy spectrum of each flavor of neutrinos is expected to be 
thermal near the surface of decoupling from matter, but then it changes due 
to propagation effects. One of these effects is scattering. Numerical modeling 
indicates that, after scattering right outside the decoupling region, neutrinos 
of a given flavor w (w = e, /i, r) have energy spectrum [34]: 



where E is the neutrino energy, L w is the (time-integrated) luminosity in the 
species w and E Qw is the average energy of the spectrum. The quantity a w is 
a numerical parameter, a w ~ 2 — 5 [34]. The non-electron neutrino flavors, z/ M , 
u T , and v T (each of them denoted as v x from here on), interact with matter 
more weakly than v e and u e , and therefore decouple from matter in a hotter 
region. This implies that at decoupling u x has a harder spectrum: Eq x ^ Eq e . 
Numerical calculations confirm this, but still leave open the question of how 
strong the inequality of energies is and of how energetic the neutrino spectra 
are. The data from SN1987A are not conclusive on this, as it will appear later. 
Indicative values of the average energies are: E QS ~ 12—18 MeV, E Qx ~ 15 — 22 
MeV. Here I consider antineutrinos only, since the v e species dominates a 
detected signal in water, and the contribution of other species is negligible for 
both the SN1987A data and for a detection of the DSNz/F. 

The second important effect of propagation on neutrinos is that of flavor con- 
version (oscillations). Conversion occurs at matter density of ~ 10 3 g • cm -3 
or smaller (see e.g. [35]), where scattering is negligible, due to the interplay 
of neutrino masses, flavor mixing and coherent interaction of neutrinos with 
the medium [36] Thus, the flux of v e of energy £ in a detector is a linear 
combination of the original fluxes in the three flavors: 



where I take into account the redshift z: here E' = E(l + z). The factor 
Pu, B {E, z) is the probability that an antineutrino produced as v w is detected as 
v e ] it describes the conversion inside the star and in the Earth and depends on 
the neutrino mixing matrix and mass spectrum. In particular, the conversion 






(3) 
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inside the star depends on the mass hierarchy (i.e. the sign of the atmospheric 
mass splitting, Amf-J and on the mixing angle 13 (assuming the standard pa- 
rameterization of the mixing matrix, see e.g. [37]). For inverted mass hierarchy 
(Am^ < 0) the propagation of antineutrinos is adiabatic if sin #13 ^10 
[35,38], resulting in a complete permutation of fluxes: Pp S + P fS = 1. In all 
the other cases (normal mass hierarchy and/or smaller 13 ) the permutation 
is only partial. I refer to the literature for more details [35,38,39]. 

The contribution of individual supernovae at different redshifts z to the DSNz/F 
at Earth is determined by the cosmic rate of supernovae Rsn{z), defined as 
the number of supernovae in the unit of (comoving) volume in the unit time. 
The rate at present is Rsn{Q) ~ O(10~ A ) Mpc -3 yr _1 . Observations as well 
as theory [40] indicate that this rate increases with the redshift, meaning that 
supernovae were more numerous in the past (Sec. 4). 

By combining the flux from an individual supernova with the rate of super- 
novae one finds the flux of v e (differential in energy, surface and time) in a 
detector at Earth: 



(see e.g. [30]). Here Q m and Q\ are the fraction of the cosmic energy density 
in matter and dark energy respectively; c is the speed of light and Hq is the 
Hubble constant. 

The expression (4) is an approximation, because it does not include a number 
of potentially relevant - but not well known - effects. One of these is the 
individual variations in the neutrino fluxes emitted, from one star to another, 
due to different progenitor mass and type, different amount of rotation and of 
convection, etc.. How large these variations can be is still an open question. 
Experimentally, an answer will come by comparing the SN1987A data with 
those of a future nearby supernova. Theoretical studies are not systematic 
enough to give a correction to Eq. (4). They hint toward little variation in the 
spectral shapes of the neutrinos, with possible variations by up to a factor of 
two in luminosity depending on the mass of the progenitor star [41]. A second 
effect not included in Eq. (4) is the deviation from the continuum limit in the 
supernova rate, due to supernova explosions in a radius of few Megaparsecs 
from Earth. These could influence the number of events recorded in the space 
of a few years [2] 3 and constitute an indirect motivation for studying the flux 
in the continuum limit: indeed, to be able to subtract the contribution of this 
flux from an observed signal would be important to conclude about a possible 
supernova event in our galactic neighborhood. 



3 I am grateful to A. Gruzinov for directing my attention to this. 
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3 Neutrino flux from an individual supernova: SN1987A 

As follows from Sec. 2, the DSNz/F in a detector depends on two sets of 
parameters. The first refers to the original spectra of the neutrinos and to 
conversion effects, while the second describes the SNR function, Rsn{z)- In 
this section I discuss the constraints on the first set of variables by analyzing 
the SN1987A data. These constraints will be used later in the calculation of 
the DSNz/F . 



3.1 The data analysis 

As input, I adopted the twelve data points from Kamiokande-II [3,4] and the 
eight events from 1MB [5], with their errors as published. I assumed that 
all these events are due to the inverse beta decay v e +p — > n + e + . The 
distributions of the observed positron energies at the two detectors are given 
in Fig. 1, together with the energies and errors of the individual events. I took 
the distance to SN1987A to be d$ 7 = 50 kpc, and neglected the ~ 10% error on 
it [42] . The effect of this uncertainty on the DSNz/F is negligible compared to 
the much larger errors of other origin. The procedure to calculate the signal at 
Kamiokande-II and 1MB given a set of parameters follows that of ref. [43]. The 
parameters subject to scan were five: two luminosities, L s , L x , two average 
energies E g, E 0x , and the mixing angle #13. I have assumed tan 2 9 12 = 0.45, 
Am^ = 8 ■ 10~ 5 eV 2 , as it is given by solar neutrinos and KamLAND (see 
e.g. [44,45]), lAmfj-J = 2.5 • 10~ 3 eV 2 from atmospheric neutrinos and K2K 
[46,47,48], and the inverted mass hierarchy. The choice to fix the inverted 
hierarchy and scan over 6*13 includes effectively the case of normal hierarchy, 
since for the latter the conversion of antineutrinos is identical to that with 
inverted hierarchy and complete adiabaticity breaking in the higher density 
MSW resonance. I followed refs. [35,38] for this and other aspects of neutrino 
conversion in the star and in the Earth. For definiteness, I fixed a s = a x = 
2.3, which give a spectral shape close to Fermi-Dirac, and therefore allow a 
meaningful comparison with several other SN1987A data analyses. The final 
results for the DSNz/F remain unchanged with the change of a s and a x . 

I parameterized the density profile of the progenitor star (necessary to calcu- 
late the matter-driven flavor conversion) as p(r) = 4-10 13 (10 Km/r) 3 g • cm -3 , 
with r being the radial distance from the center. For simplicity, I ignored the 
uncertainty on p(r), as well as the errors on the values of the solar and at- 
mospheric oscillation parameters. The inclusion of these uncertainties has no 
appreciable effect on the final results for the DSNz/F. I have not included late 
time shockwave effects [49], as their impact on a time integrated signal is 
negligible with respect to the large statistical errors and other uncertainties 
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Fig. 1. The energy spectra of the positrons observed at K2 and 1MB. The energies 
of the individual events with their errors are shown as well, as in ref. [56] . 



(see e.g. [28] for this particular aspect). The experimental parameters, such 
as efficiency curves and energy resolution functions, are as in [43], and the 
detection cross section was taken from [50] (Eq. (25) there). 

Given the sparseness of the data, the maximum likelihood method of analysis is 
the most appropriate. Following Jegerlehner, Neubig and Raffelt [51], I obtain 
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the likelihood function, C 87 (E 0S , E 0x , L s , L x , sin 2 6*13), and the quantity Xs7 : 

xl 7 = -2\nC 87 . (5) 



Once the minimum of x| 7 , X%7mini nas been found, the five- dimensional region 
of parameters that are allowed at a given confidence level (C.L.) is given by: 

Axg 7 = xl 7 ( E os, E 0xj L B , L x , sin 2 13 ) - xl 7 ,min < Xs , (6) 

where Xs = 5.86, 9.24, 15.09 for 68, 90, 99% C.L.. Like in [51], here the value of 
Xg7 will be given up to a constant, which is irrelevant for parameter estimation, 
according to Eq. (6). I emphasize that the maximum likelihood method does 
not involve any binning of the data [51]. Thus, the bins used in figures 1 and 
2 are only for illustration and do not influence the likelihood function. 

The scan was performed in the five-dimensional box given by E os = 3 — 30 
MeV, E 0x = 3 - 30 MeV, L- e = (1.5 - 45) • 10 52 ergs, L x = (1.5 - 45) • 10 52 
ergs, sin 2 613 = 10 -7 — 10~ 2 . The latter interval covers all the possibilities 
of spectrum permutation due to conversion in the star [38]. No hierarchy of 
neutrino spectra was imposed a priori. Results with certain priors will be 
discussed briefly later. 

The five-parameters analysis adopted here is statistically meaningful given 
the number (20 in total) of data points available. It is original, because it 
combines generality (the absence of theoretical priors) and great detail in the 
inclusion of neutrino conversion effects 4 . It is known that a smaller num- 
ber of parameters is enough to well reproduce the SN1987A data [56] and 
thus would probably suffice to predict the DSNz^F. Still, the five-parameters 
method has been preferred here because, aside from the specific connection 
to the DSNz/F, it serves a more general purpose: to give a state-of-art answer 
to the question of what we know on the neutrino energies and luminosities 
at the production site. This question is important in connection with other 
physics inside the supernova (neutrino transport, R-process nucleosynthesis, 
etc.). Another reason for the choice is that this method includes the case when 
the observed neutrino spectrum is strongly distorted with respect to an effec- 
tive thermal distribution. This can happen if the original v e and v x spectra 
have very different average energies and luminosities (see Sec. 3.2 and Fig. 2). 



The statistical analyses in the past literature (e.g. [51,52,53,54,55,56]) used priors. 
Most of them had only two fit parameters, with fixed flavor conversion pattern and 
fixed ratios of average energy and of luminosities in v e and v x . 



8 



3.2 Analyzing Kamiokande and 1MB individually 

It is useful to consider the K2 and 1MB data sets separately first. Their energy 
spectra are shown in Fig. 2, and their average energies are summarized in Table 
1. One can see that the two spectra differ substantially. The differences are 
in good part due to different experimental settings (higher energy threshold 
for 1MB, different volumes, energy resolutions and detection efficiencies). Still, 
on top of these technical differences, a tension exists between the two spectra 
5 . Specifically, the fact that the 1MB spectrum has a maximum at about 40 
MeV contrasts with the K2 signal, which is compatible with an exponentially 
decreasing spectrum in that energy region. 

In the light of this tension, it is meaningful to study the x 2 functions of K2, 
x\ 2 i an d 1MB, x 2 imbi separately to understand how each data set influences 
the combined function Xs7- 

As expected given the large number of fit parameters, both x\i an d x\mb have 
many degenerate minima, that correspond to the same spectrum of events. 
Two of these minima are given in Table f for illustration. From these sets of 
minima, one can see that: 

• The K2 observed spectrum is best reproduced by a superposition of two 
original spectra, one very soft, Eq ~ 4—5 MeV and the other more energetic: 
E ~ 12 — 14 MeV. The soft component has to have very high luminosity, 
L ~ 3- 10 53 ergs, higher by a factor of several with respect to the hard one. In 
this composite spectrum, the soft part accounts for the peak of the K2 data 
in the lowest energy bin, E ~ 10 MeV, while the harder part reproduces 
the tail of the observed spectrum at higher energy. This is illustrated in 
Fig. 2. The scenario favored by the K2 data can be realized for normal 
hierarchy, or inverted hierarchy with small #13, a soft v e original spectrum 
and a harder u x one (the case with soft v x and hard v e fits equally well, 
but is not theoretically motivated). One must be aware, however, that this 
spectrum contrasts with the theory (Sec. 2) because the soft component 
has too low energy and too high luminosity, comparable to the total energy 
output predicted for a supernova 6 . 



5 The angular distribution of the events at both detectors is only marginally com- 
patible with the predictions for inverse beta decay data, suggesting the presence of 
an anomaly in the data [57]. The standard interpretation of the SN1987A signal in 
terms of inverse beta dacay is not excluded, however, and therefore it is adopted 
here. 

6 The conclusion of ref. [43], that the K2 data alone do not favor a composite spec- 
trum, is not in conflict with the results of this work. It reflects the more conservative 
assumptions adopted in [43] about the original neutrino fluxes, such as the condition 
of comparable luminosities: L x /L s = 0.5 — 2. 
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Table 1 



Summary of relevant quantities in the points of maximum likelihood of K2 and 
1MB individually, as well as in the point of maximum combined likelihood and in 
the 68% C.L. allowed region. The three specific points used here were chosen for 
illustration among the degenerate maxima of the likelihood function (minima of x 2 )- 
The values of x 2 are given up to a constant (irrelevant for parameter estimation). 
Where meaningful, the values given by the data are presented too (from ref. [43]). 



• The 1MB data favor a purely thermal spectrum - instead than a superpo- 
sition of thermal spectra - with average energy E ~ 13 — 15 MeV and 
luminosity L ~ 0.5 ■ 10 53 ergs, in acceptable agreement with theory. These 
parameters are determined by the average energy of the data and by the 
number of events. A thermal spectrum is favored over a composite one be- 
cause of its smaller width (see e.g. [43] for a detailed discussion), that better 
reproduces the narrow spectrum of the 1MB data. The spectral character- 
istics favored by 1MB are realized if the mass hierarchy is inverted with 
adiabatic conversion (large #13), for which the whole v e signal is due to the 
original u x flux. The other cases, where the conversion is only partial, fit 
equally well if the original v x and v e fluxes have the same spectrum or if 
one of them is suppressed with respect to the other one. The suppression 
can be due either to small luminosity or to average energy much below the 
1MB detection threshold [Eth ~ 20 MeV [5]). The latter point is especially 
important, as will be seen. 
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Fig. 2. The observed energy spectra of events at K2 and 1MB, as in Fig. 1, compared 
with the predicted spectra in the points of minimum x 2 f° r K2 only, 1MB only and 
combined K2 and 1MB data sets. The values of the parameters in these points are 
given in Table 1. 

3.3 Combined analysis: results 



Comparing the neutrino spectra favored by K2 and by 1MB separately, one 
infers that a good combined fit exists. The key to see this is to observe that 
the average energy favored by 1MB is similar to that of the hard component 
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of the spectrum favored by K2 (see Table 1). Moreover, a very soft spectral 
component would affect 1MB only marginally, due to threshold effects, as was 
mentioned in Sec. 3.2. Thus, one expects that the combination of the two data 
sets will favor a composite soft+hard spectrum qualitatively similar to that 
favored by the K2 data. The closeness to the K2-only result is motivated also 
by the fact that the K2 data dominate the statistics. 

Some features of the xh function can be predicted as well. One expects de- 
generacies, even though not as extended as in the case of K2 and 1MB sep- 
arately. In particular, a degeneracy in \\j should exists between the region 
with Eos < Eq x and that with Eq s > Eq x . The reason lies in the fact that, if 
the neutrino conversion is only partial (probability < P xS < 1), the neutrino 
spectrum at Earth is symmetric under the transformation P ss dN s (E')/dE' <-> 
(1 — P ss )dN x (E') / dE' (see Eq. (3)). Such symmetry implies an approximate 
symmetry (broken by effects of oscillations in the Earth and by the energy 
dependence of P g g) between the E os < E 0x and E os < E Qx portions of the 
parameter space. 

The likelihood analysis, summarized in Table 1 and Figs. 2 and 3, confirms 
the intuitions. The predicted degeneracy is observed. The Xs7 nas two minima 
in the points (log(sin 2 13 ), E os , E 0x , L s , L x ) = (-7,4.6 MeV, 12.7 MeV, 3.4 • 
10 53 ergs, 0.51 • 10 53 ergs) ( X 2 87 = 84.1) and (log(sin 2 13 ), E os , E 0x , L g , L x ) = 
(-3.6, 12.5 MeV, 5.2 MeV,6-10 53 ergs, 1.7-10 53 ergs) (xl 7 = 84.7). They corre- 
spond to the type of spectrum preferred by K2. Here, the same considerations 
done for the K2-only fit apply regards the conflict between the best fit scenario 
and theory. It is important to stress, however, that the Xs7 function is very 
shallow, with many points having essentially the same goodness of fit as the 
minimum. The allowed region in the parameter space is very extended and 
includes theoretically motivated scenarios at 68% C.L. (Fig. 3). The region is 
compatible with the constraints on the neutrino spectrum put by light element 
synthesis during galactic chemical evolution [58,59] (see also the related works 
[60,61]), which favor rather soft muon and tau neutrino spectra: Eq x ~ 15 — 18 
MeV [58]. A detailed comparison of those constraints with the results of this 
paper is not possible, however, due to the different working assumptions used 
here. 

The upper panels of Fig. 3 illustrate the details of the allowed region. There 
I give the projections on the Eq s — Eq x and Eq x — L x planes of the 5D al- 
lowed regions given by Eq. (6) for different confidence levels, together with 
the projections of the points of maximum likelihood. The figure confirms that 
a rather large portion of the parameter space is allowed, with better sensitiv- 
ity to the v x flux with respect to the v e one. Unless assumptions are made 
on the oscillation pattern, the v e flux is unconstrained, due to the possibil- 
ity that the mass hierarchy be inverted with perfectly adiabatic conversion 
(sin 2 6 ^ few ■ 10~ 4 ) inside the star. Under such conditions the original u e 
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Fig. 3. Projections of the 68%,90%,99% C.L. regions allowed by the SN1987A data 
on the planes Eq s — Eq x and Eq x — L x , without any prior (upper panels) and with 
the hierarchy Eq s < Eq x (lower panels). The dots in each panel mark the projections 
of the points of maximum likelihood (see Table 1). The entire plane Eq s — L s (not 
shown) is allowed at 68% C.L.. 

flux is completely converted into v x and does not affect the observed v e signal. 
The v x flux is constrained loosely: every value of E 0x gives a good fit provided 
that the other parameters are suitably adjusted, and a similar statement is 
valid for L x . Notice that the v x luminosity and average energy are not con- 
strained from below: indeed, they are allowed to be zero, in the case when 
the v e spectrum alone gives a good fit to the data. This can happen for all 
conversion patterns except the one with complete v e — v x permutation, and 
requires E os ~ 13 — 15 MeV. In all cases, the region with E os ^ 5 MeV and 
E 0x ^ 5 MeV is excluded because it gives a very soft neutrino spectrum, with 
a too small (or even vanishing) number of events at 1MB. 

If the analysis is restricted to the region motivated by the theory, E os < E 0x 
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Fig. 4. The best fit point and the regions allowed by the SN1987A data at 
68%, 90%, 99% C.L. in the space of the average energy and luminosity of the v e 
flux at Earth (after conversion in the star). For comparison, I show the 99% C.L. 
region found in ref. [56] (dashed contour). There the same flux was described by a 
spectrum of the form (2). 

(lower panels of Fig. 3), some conclusions change: Eq x must be larger than 
~ 10 MeV and E os can not exceed ~ 16 MeV. 

The Table 1 gives further details on one of the best fit points, confirming 
that it reproduces the data well in the number of events, Nk2 and Njmb and 
average energies, {E) K2 and (E) IMB . The goodness of the fit also appears 
from the superposition of the best fit spectrum with the data, in Fig. 2. 

It is informative to present the results in a complementary way: by using a 
smaller number of variables and choosing parameters that describe the v e flux 
at Earth, instead than at the production point. Here I choose the average 
energy and luminosity of such flux, E av and L tot . Fig. 4 shows the allowed 
region in the space of E av and L to t, obtained after marginalizing the likelihood 
function over the three parameters orthogonal to these. For comparison, in the 
figure I also show an example of result obtained by fitting the data with two 
parameters only (i.e., keeping the other three parameters fixed instead than 
marginalizing over them). The example is taken from ref. [56] (fig. 2 there), 
where the v e flux at Earth is described by a spectrum of the form (2). The 
marginalized region confirms the finding of a particularly soft and luminous 
P e flux being favored by the data. The region is compatible with but more 
extended than the two-parameters result from ref. [56]. This reflects the fact 
that a five-variables parameterization includes neutrino spectra that fit the 
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data well but are not reproduced with a smaller number of parameters. It is 
confirmed, therefore, that five-dimemsional and the two-dimensional methods 
are not equivalent. 



4 The supernova rate 

The cosmic rate of supernovae, Rsn{z), can be inferred from astrophysical 
data in a variety of ways. Here I review two: 

(1) Measurements of the SNR from observations of core collapse supernovae. 
This is the most direct method. To date, four measurements of the SNR 
have been made in this way, covering the interval of redshift z ~ — 0.9 
[14,15,16]. They are summarized in Table 2 and in Fig. 6. While the con- 
nection between the observations and the SNR is immediate, in principle, 
one should keep in mind that the results are affected by dust obscura- 
tion (extinction) and possible misidentification of the observed objects. 
These effects can be modeled theoretically and subtracted. Three of the 
measurements in Table 2 have been corrected in this way. The point at 
z = 0.26, indicated with an empty circle in Fig. 6, is not corrected for 
extinction and for this reason it will not be included in the analysis here. 

(2) measurements of the star formation rate (SFR), R SF (z). From these, the 
SNR is found through the equation: 



(see e.g. [30]) where <f)(m) is the initial mass function, decreasing roughly 
as a power —2 of the mass m, and M is the mass of the Sun, M ~ 1.99 • 
10 30 Kg. The limits of integration represent the interval of mass for which 
a star becomes a core collapse supernova. With respect to direct searches 
of supernovae, this method has the advantage of higher statistics, since 
several measurements of the SFR are available, extending up to z ~ 6 — 7 
(see e.g. [17,18,19] and [62] for an overview and further references). On the 
other hand, the SNR obtained with this approach is indirect: it is affected 
by a number of theoretical uncertainties through the relation (7), such 
as those on the initial mass function and on the interval of progenitor 
masses. The lower mass cut in (7) influences the normalization of the 
SNR strongly. I refer to [62] for an in-depth discussion of these aspects. On 
top of the uncertainties mentioned, the SFR measurements are affected 
by extinction, like the direct supernova observations. Likely, the largest 
theoretical error is associated to the normalization of the SNR, however 
the possibility that the ratio of SFR and SNR be redshift-dependent - 
thus violating Eq. (7)- is not excluded. This is a possible source of error 
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0.59 ± 0.24 


[16] 


0.26 




[15] 


0.3 (average) 


o ci +0.88 
Z - Oi -0.75 


[15] 


0.7 (average) 


Q OR+1.03 
°- yo -1.06 



Table 2 

The measured supernova rate at different redshifts. The errors represent the 68% 
confidence level intervals. The values of the redshift in the third and fourth row are 
averages over bins of redshift. The point at z = 0.26 is not corrected for exctinction, 
while the other three points are. 

as well. 

Clearly, the two methods are complementary. Their results agree in the basic 
facts: the SNR today is of order Rsn{Q) ~ 10 -4 yr _1 Mpc~ 3 . The rate increases 
with z and is consistent with a broken power law: 



R SN (z) = R SN (0)(l + z)P for z<\ 

= R SN (0)2 p ~ a (1 + z) a for z > 1 . (8) 

The values of Rsn{0) favoured by the two approaches are somehow different, 
with method (2) giving a rate larger by a factor of ~ 2 — 3 with respect to 
method (1). This is probably a manifestation of the uncertainties that affect 
both methods. For example, it was checked that the discrepancy disappears if 
the lower mass cut in (7) is increased to 1OM [62]. 

In this work I use method (1), motivated by its being direct and therefore 
more robust, and suitable for a conservative estimate of the DSNVF. For this 
specific application, its being limited to small resdhift is of little consequence, 
since the DSNVF above realistic energy thresholds is dominated by super- 
novae closer than z ~ 1 [30]. For this reason, here I fix a = 0, a value favored 
by the SFR data [62]. I perform a maximum likelihood analysis of the three 
extinction-corrected measurements of the SNR given in Table 2, to find the 
allowed region of the parameters Rsn{®) and (3 of Eq. (8). In the approxima- 
tion of gaussian errors, I calculate the likelihood function £snr(Rsn(0), P) 
and Xsnr = -2\n£ S Nn, 

The results are shown in figs. 5 and 6. The likelihood is maximal in the 
point R S n(0) = 0.67 • 10~ 4 yr^Mpc" 3 , (3 = 3.44, where the value of Xsnr is 
XsNRmin = 3.59. Fig. 5 shows the point of maximum likelihood and the con- 
tours defined by Xsnr. ~~ XsNRmin = 2-3, 4.61, 6.17, corresponding to 68.3, 90, 
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95.4% confidence level (C.L.) 7 . The results of the likelihood analysis confirm 
that the fit of the data with a power law curve is satisfactory. The allowed 
region is rather extended, with R~sn{Q) and (3 varying by a factor of six and 
three respectively. In Fig. 6 the SNR functions allowed at given confidence 
level (shaded bands) and the best fit function are compared with the input 
measurements of Table 2. The figure also gives the SNR function found in ref. 
[62] using method (2), with the modified Salpeter B initial mass function (Ta- 
ble 2 in [62]) and the interval of progenitor masses as in Eq. (7). Confirming 
what already mentioned above, this curve has very similar logarithmic slope, 
(3 = 3.35, but normalization higher by a factor of ~ 2 with respect to the best 
fit function found in this work. 
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Fig. 5. Best fit point and isocontours of x 2 i n the space of the parameters de- 
scribing the SNR function, Rsn(z). These are the intercept, Rsn{0) (in units of 
10" 4 yr" 1 Mpc" 3 ) and the power, 0. The contours refer to 68.3, 90, 95.4% C.L.. 

In connection to this, one may consider performing a global analysis of the 
direct SNR measurements and of the SFR data together. Such analysis war- 
rants a separate work, necessary to make the different data sets compatible 

One should keep in mind that these confidence levels are not exact, as they apply 
rigorously only to likelihood functions that are perfectly gaussian in the parameters 
of interest. 
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Fig. 6. The SNR (in units of 10 -4 yr _1 Mpc -3 ) as a function of the redshift, z. 
The colored (shaded) bands are the families of curves having parameters within the 
68.3, 90, 95.4% C.L. contours of Fig. 5. The lower dashed line refers to the best 
fit parameters of my analysis. The measurements of the SNR, with 1 a error bar, 
are shown (from Table 2). Of these, the point marked with an empty circle was 
not included in the analysis because it is not corrected for exctinction. The upper 
dashed line results from using, instead than core collapse data, the SFR data with 
a specific initial mass function (from ref. [62], see text). 

with each other; this in turn requires, among other things, a careful study 
of the connection between SNR and SFR, and an evaluation of the different 
methods used by different authors in their statistical analyses and in the cor- 
rection for dust obscuration. This has not been done so far and represents 
the next step with respect to this paper. In consideration of the uncertain- 
ties on the normalization of the SFR, I expect that the combination with the 
SFR data would improve the constraint on j3, with only little improvement on 
Rsn(0). 

Studies have indicated that the z = measurement used here may under- 
estimate the local SNR by a factor of two or so [63,2]. In absence of a new 
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determination of the local rate, one can only give examples of how the likeli- 
hood results would change if the z = point was higher. If I rescale this point 
and its error by a factor of two (three) the new point of maximum likelihood is 
R SN (0) = 1.26 • 10~ 4 yr^Mpc" 3 , f3 = 2.21 (R SN (0) = 1-72 • 1(T 4 yr^Mpc" 3 , 
(3 = 1.55), with only minor change in the value of XsNRmin- As wm be shown, 
the effect of a higher Rsn{Q) is partially compensated by that of a smaller (3, 
so that the impact of this change on the DSNVF is only of tens of per cent. 



5 Diffuse neutrino flux 

5.1 The calculation 

In the framework adopted here, the DSNz/F depends on seven parameters: five 
describing the neutrino flux (two average energies, two luminosities and one 
mixing angle) and two describing the SNR (the rate at z = and the power, 
(3). To calculate the diffuse flux and the uncertainty on it, it is necessary 
to combine the information from both the SN1987A and the SNR data sets 
consistently. 

I perform this combination as follows: 

• I obtain the total likelihood of the neutrino and SNR data together, £ tot . 
Since the two sets of data are uncorrelated, this is done by simply multiply- 
ing the two individual likelihoods. The combined x 2 is the sum of the two 
pieces: x 2 = xli + Xsnr- 

• Using Eq. (4) I compute the likelihood for the v e diffuse flux in a detector, 
$, by marginalizing C tot with respect to all the six quantities orthogonal to 
$. If the likelihood function is close to a gaussian near its maximum, for a 
given $ the marginalized x 2 ~ called x% from here on - is well approximated 
by the global \ 2 minimized with respect to the other six variables. Here I 
adopt this approximation, widely used in data analysis literature (see [64] 
for an example). 

• I use x% t° fi n d the best value of $, defined as the one that minimizes x%i 
and the 99% C.L. interval for it, defined by the difference x%~ X% mm = 6.65. 

Results are given for three cuts in the neutrino energy: E > 19.3,11.3,5.3 
MeV. The first corresponds to the threshold applied in the search of the 
DSN^F at SK [6], while the second represents a potential improvement that 
SK could achieve with the proposed Gadolinium addition [65]. The last cut, 
E > 5.3 MeV, is the technical limit given by the photomultipliers distribu- 
tion in the SK tank and would apply only in the - currently unfeasible - 
case of complete background subtraction. It should be stressed that with such 
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low threshold the DSNi/F starts to have a non-negligible contribution (about 
~ 50% [30]) from supernovae at z > 1, and thus it depends on the value of 
the power a in Eq. (8). Here I used a = 0, therefore the results of this work 
for the lowest threshold have indicative character only; they could change by 
several tens of per cent for different values of a. 

The effects of neutrino flavor conversion inside the star are included. I calcu- 
lated the effects of oscillations inside the Earth, using a realistic matter density 
profile [66], and in the assumption of isotropic DSN^F. I find that these os- 
cillations affect the DSN^F by less than ~ 10%, and therefore I neglect them 
for simplicity. 

The same procedure of marginalization discussed for the flux can be done for 
the number of events in a detector. Here I do the calculation for SK and the 
inverse beta decay events: v e +p — > n + e + . This is by far the dominant signal 
in water, which justifies neglecting other reactions. I took the fiducial volume 
of SK of 22.5 kilotons and detection efficiency of ~ 93% above ~ 7 MeV [3,51]. 
In reality, the efficiency depends on the specific experimental cuts. In the case 
of the published SK analysis the efficiency is 47% (79%) below (above) 34 MeV 
[67]. Thus, considering the exponential decay of the flux with the energy, the 
rates given here would have to be rescaled down by roughly a factor of 2 to 
be applicable to the current SK setup. 

5. 2 Results 

The results for the flux and numbers of events are presented in Table 3 and 
Fig. 7. 

For the current SK threshold of 19.3 MeV, the v e flux in a detector is of the or- 
der of ~ 10" 1 cm~ 2 s _1 . The value of maximum likelihood is $ = 0.15 cm~ 2 s~ 1 
and is obtained with the parameters that minimize both Xs7 an d Xsnri Sees. 
3 and 4. With 99% C.L., the flux must be larger than ~ 0.05 cm~ 2 s _1 and can 
not exceed 0.35 cm _2 s _1 , almost a factor of 4 below the SK limit. The event 
rate is below 0.7 events/year. The width of the 99% interval is about a factor 
of 7-8 both in the flux and event rate. 

The flux increases by almost one order of magnitude if the threshold is lowered 
to 11.3 MeV. This is explained by the exponential decay of the flux with 
energy. The event rate instead increases more moderately with the lowering 
of the threshold. This is because the detection cross section is proportional to 
the square of the energy, a oc E 2 , and therefore the spectrum of the observed 
positrons decays less rapidly than the neutrino spectrum. From Table 3 it 
appears that for threshold at 11.3 MeV or lower the goal of one event/year is 
possible, but not guaranteed. 
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Table 3 

The predicted flux of u e in a detector in the point of maximum likelihood and in 
the intervals of 68, 90, 99% C.L. for different energy thresholds. The numbers in 
brackets have the same meaning, but for the rates of v e +p — > n + e + events at 
SuperKamiokande for ~ 93% detection efficiency (see text). 

To expect at least one event/year with 99% C.L. a larger detector than SK is 
necessary. A volume of water 2.5 times the volume of SK would be sufficient 
for the lowest threshold of 5.3 MeV, while with the current SK threshold a 
detector ~ 10 times larger is required. This may become a reality with the 
next generation Megaton detectors. With a fiducial volume 20 times larger 
than SK, these detectors would register ~ 2 — 14 events/year for the higher 
threshold, with a maximum of ~ 44 events/year with detection above 5.3 
MeV. 

What is the contribution of the individual sources of error to the total un- 
certainty on the DSNz/For event rate? Fig. 7 answers this question for the 
event rate. For the three thresholds, the figure shows the error bars obtained 
by marginalizing the likelihood over a subset of parameters while keeping the 
others fixed at their best fit values. It appears that the error bar from the SNR 
(with fixed neutrino parameters) is smaller than the one due to the neutrino 
spectra (e.g. fixed SNR). The first is about a factor of 3 wide, while for the 
second the width is a factor of 4-7 depending on the threshold. Notice that 
this width is larger for higher threshold, reflecting the larger uncertainty in 
the high energy tail of the SN1987A data with respect to their full spectrum. 

Let me now discuss the validity of the results and possible generalizations. 
A relevant question is how the interval of values of the DSNi/F depends on 
the details of the analysis of the SN1987A data, and, in particular, on pos- 
sible priors on the neutrino parameters. I checked that the results have lit- 
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Fig. 7. The number of P e +p — > n + e + events per year at SuperKamiokande 
as a function of the energy threshold Eth (horizontal axis). Results are given for 
E t h = 5.3, 11.3, 19.3 MeV in neutrino energy. For each of these, three cases are 
shown (shifted in energy for visual convenience): (i) total error, from the marginal- 
ization of the likelihood function in the whole parameter space, (ii) error from the 
neutrino spectrum only, obtained by marginalization with the supernova rate func- 
tion fixed at its best fit point, (iii) error from the supernova rate only, obtained by 
keeping the neutrino parameters fixed at the best fit values. The markers give the 
values in the best fit point, while the bars represent the 99% C.L. interval ("error"). 

tie dependence on those. For example, the interval for the flux above 19.3 
MeV (see Table 3) is unchanged if I impose (simultaneously) the constraints: 
10 MeV < E os < 20 MeV, E os < E x and 0.5 < L s /L x < 2. This is explained 
with the parameter space being largely degenerate, so that a restricted portion 
of it still covers almost all the physically different possibilities. For the same 
reasons the results are also practically insensitive to different parameteriza- 
tions of the original fluxes (e.g., different a s and/or a x ) 8 . 

The analysis here does not include systematic errors. Still, here I briefly discuss 
how the prediction of the DSNVF would change if some of the inputs turned 
out to be systematically wrong and were corrected. Let us consider the event in 
which i?57v(0) is higher than what used here (see discussion in Sec. 4). With 
i?5^v(0) increased by a factor of two (three) the maximum likelihood value 
of the flux is $ = 0.20 cm~ 2 s" 1 ($ = 0.23 cm~ 2 s~ 1 ) above the current SK 
threshold of 19.3 MeV. This is at most a ~ 50% increase with respect to the 
result in Table 3. Another example refers to the poorly known normalization 
of the SNR. Specifically, my results could become a factor of two higher if it 
is confirmed that normalization of the SNR is two times higher as favored by 
the SFR measurements (Sec. 4 and Fig. 6). A combined fit of the supernova 

8 Clearly, different parameterizations of the original fluxes give differences in the 
region of the energies and luminosities allowed by SN1987A, however. 
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data and SFR measurements would reduce the uncertainty on the slope (3 and 
therefore the total error on the DSNz/F. To check how this error would change, 
I performed the same calculation outlined in this section with j3 fixed at its 
best fit value, f3 = 3.44, which happens to be nearly the same both in the SNR 
and in the SFR data analyses (Sec. 4). I find that the results for the DSNz/F 
do not change appreciably, due to the large degeneracies between different sets 
of parameters. 

Besides technical details, it has to be stressed that the results given here are 
valid if the SN1987A neutrino flux is typical and thus represents a generic 
output of a core collapse supernova. This question will be answered by data 
from a future galactic supernova. At the moment, recent calculations show that 
the main observational features of SN1987A are reproduced by the standard 
neutrino-driven explosion mechanism, disfavoring possible anomalies in the 
SN1987A event [68]. 



6 Discussion and conclusions 

Let me summarize this work. I have calculated the v e component of the DSNz/F 
in a detector, using the information on the neutrino spectra from SN1987A and 
the information on the supernova rate from direct observations of core collapse 
supernovae. I calculated the likelihood functions for the SN1987A data and for 
the supernova rate measurements, and marginalized the combined likelihood 
to find the interval of DSNz/F allowed at a given confidence level. 

The SN1987A data favor a composite neutrino spectrum, meaning that a 
scenario with inverted mass hierarchy and sin 2 # 13 ^ 10~ 4 is disfavored (in 
agreement with refs. [52,53]). The best fit spectrum reproduces both the K2 
and 1MB data well, even if 1MB alone favors a thermal spectrum over a com- 
posite one. It is characterized by a very soft and very luminous component, 
with parameters in contrast with general theoretical arguments and with nu- 
merical simulations. The allowed region of the parameter space, however, is 
very extended and includes more natural scenarios at 68% C.L.. 

The supernova rate measurements allow a present rate of about (0.2 — 1.2) • 
10 -4 Mpc~ 3 yr" 1 and a power of increase with z of about 2 - 6 . Also in this 
case a large region of parameters is allowed. The function that best fits the 
data has similar power, (3 ~ 3, but different normalization (about a factor of 
two smaller) with respect to what is inferred using data on the star formation 
rate and the proportionality of this rate to the supernova rate. 

The results for the DSNz/F {p e component) show that the flux above the cur- 
rent SK threshold is likely to be about one order of magnitude below the 



23 



current upper limit of 1.2 cm _2 s _1 , and is smaller than this limit by a factor 
of ~ 4 with 99% confidence. This factor of 4 is the minimal improvement 
that SK should achieve if the energy threshold remains the same. Any lower- 
ing of the threshold would significantly enhance the possibility to detect the 
DSNz/F. Still, however, if my prediction is correct, a Megaton detector would 
be necessary to expect a detection with very high confidence. 
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Fig. 8. The results of different authors (from left to right: refs. [25], [26] (upper 
bound only), [31,30], [29] and this work) for the flux of v e above a threshold of 19.3 
MeV (in neutrino energy). For this work the colored (shaded) bands correspond to 
68%, 90% and 99% confidence levels, with the central like marking the maximum 
likelihood value. The SK limit is shown as well (dashed line). 

How does this work compare to others in the field? Fig. 8 answers this question, 
by presenting the results of different papers [25,26,31,30,29] for the v e diffuse 
flux above 19.3 MeV, with the interval of uncertainty associated to it. The 
authors quoted in the figure are only a small sample of all those that have 
calculated the DSNz/F; they were chosen for illustration purpose and because 
their results can be directly compared to each other's and to mine. Each 
of them used neutrino spectra from a number of numerical simulations, and 
inferred the supernova rate from the star formation rate. They also gave a 
tentative uncertainty on their result, associated with the lack of consensus 
on numerical codes and with the error on the star formation rate. The error 
bars in the figure, as well as the central points, do not have a statistical 
meaning except for the result of this paper, where the shaded areas represent 
the intervals allowed at 68%, 90% and 99% CL. and the central line marks 
the maximum likelihood value. 

From the figure it appears that the interval of flux calculated here overlaps 
with the prediction of other authors within the error, but extends to lower 
values. In particular, it is the only calculation which allows a flux smaller than 
~ 0.1 cm~ 2 s _1 . The differences between my results and the larger flux - close 
to the SK limit - allowed by Strigari et al. [29] and by Ando and Sato [31,30], 
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are due to the different inputs used by these authors: a factor of about two 
higher SNR (see discussion in Sec. 4) and larger neutrino flux at high energy, 
due to the harder neutrino spectra used, with E 0x up to ~ 24 MeV. These 
spectra are motivated by numerical calculations of neutrino transport, but are 
known to be at best marginally compatible with the SN1987A data [51,53,54], 
a fact that this paper confirms 9 . Hartmann and Woosley [25] used a supernova 
rate function similar to that of Strigari et al. and of Ando and Sato, together 
with softer neutrino spectrum: a thermal spectrum with temperature of ~ 4 
MeV. This explains their finding a lower flux with respect to other authors. 
The result by Kaplinghat et al. [26] is a conservative upper bound obtained 
with a unoscillated neutrino spectrum and the maximum SNR compatible 
with metal enrichment history. 

The comparison with the previous literature shows that this work is comple- 
mentary to others. It is especially useful because it gives an idea of how the 
predicted DSNz/F can change with a change of the method of calculation or of 
the inputs. This should be encouraging to work to improve our knowledge of 
the inputs, in particular of the neutrino spectra and of the normalization of 
the SNR. This paper also makes the point that the flux may be smaller than 
generally expected, thus supporting the case for a next generation of neutrino 
detectors with larger volumes and/or lower energy thresholds. 



Acknowledgements 

I am very grateful to A. Friedland, W. Haxton and A. Yu. Smirnov for helpful 
discussions and comments on the first version of the manuscript. I warmly 
thank E. Cappellaro, M. Giavalisco, M. Maltoni and M. Turatto for fruitful 
private communications. I am indebted to the anonymous referee for crucial 
comments. I acknowledge support from the INT-SCiDAC grant number DE- 
FC02-01ER41187. 



References 

[1] N. Arnaud et al., Astropart. Phys. 21, 201 (2004). 

[2] S. Ando, J. F. Beacom and H. Yuksel, Phys. Rev. Lett. 95, 171101 (2005). 

9 This poor or no compatibility is not surprising. Indeed, numerical calculations 
are currently limited to rather crude neutrino transport or to simulating only the 
first second or so after the core bounce. Therefore, they are not directly comparable 
to the SN1987A data. 



25 



[3] K. Hirata et al., Phys. Rev. Lett. 58, 1490 (1987). 

[4] K. S. Hirata et al., Phys. Rev. D38, 448 (1988). 

[5] R. M. Bionta et al., Phys. Rev. Lett. 58, 1494 (1987). 

[6] M. Malek et al., Phys. Rev. Lett. 90, 061101 (2003). 

[7] K. Eguchi et al., Phys. Rev. Lett. 92, 071301 (2004). 

[8] W. C. Haxton and C. W. Johnson, Nature 333, 325 (1988). 

[9] J. N. Bahcall, Nature 333, 301 (1988). 

[10] C. K. Jung (1999), hep-ex/0005046. 

[11] R. J. Wilkes (2005), hep-ex/0507097. 

[12] K. Nakamura, Int. J. Mod. Phys. A18, 4053 (2003). 

[13] L. Mosca, Nucl. Phys. Proc. Suppl. 138, 203 (2005). 

[14] E. Cappellaro et al, Astron. and Astrophys. 351, 459 (1999). 

[15] T. Dahlen et al, Astrophys. J. 613, 189 (2004). 

[16] E. Cappellaro et al., Astron. and Astrophys. 430, 83 (2005). 

[17] S. Cole et al, Mon. Not. Roy. Astron. Soc. 326, 255 (2001). 

[18] K. Glazebrook et al, Astrophys. J. 587, 55 (2003). 

[19] I. K. Baldry and K. Glazebrook, Astrophys. J. 593, 258 (2003). 

[20] D. Schiminovich et al, Astrophys. J. 619, L47 (2005). 

[21] L. E. Strigari, et al., JCAP 0504, 017 (2005). 

[22] M. Fukugita and M. Kawasaki, Mon. Not. Roy. Astron. Soc. 340, L7 (2003). 

[23] T. Totani, K. Sato and Y. Yoshii, Astrophys. J. 460, 303 (1996). 

[24] R. A. Malaney, Astropart. Phys. 7, 125 (1997). 

[25] D. H. Hartmann and S. E. Woosley, Astropart. Phys. 7, 137 (1997). 

[26] M. Kaplinghat, G. Steigman, and T. P. Walker, Phys. Rev. D62, 043001 (2000). 

[27] S. Ando, K. Sato and T. Totani, Astropart. Phys. 18 (2003) 307. 

[28] S. Ando and K. Sato, Phys. Lett. B559, 113 (2003). 

[29] L. E. Strigari, M. Kaplinghat, G. Steigman, and T. P. Walker, JCAP 0403, 007 
(2004). 

[30] S. Ando and K. Sato, New J. Phys. 6, 170 (2004). 
[31] S. Ando, Astrophys. J. 607, 20 (2004). 



26 



[32] A. G. Cocco, A. Ereditato, G. Fiorillo, G. Mangano and V. Pettorino, JCAP 
0412, 002 (2004). 

[33] J. F. Beacom and L. E. Strigari, |hep^pn"/ 0508202 

[34] M. T. Keil, G. G. Raffelt, and H.-T. Janka, Astrophys. J. 590, 971 (2003). 

[35] A. S. Dighe and A. Y. Smirnov, Phys. Rev. D62, 033007 (2000). 

[36] S. P. Mikheev and A. Y. Smirnov, Sov. Phys. JETP 64 (1986) 4 [Zh. Eksp. 
Teor. Fiz. 91 (1986) 7]. 

[37] P. I. Krastev and S. T. Petcov, Phys. Lett. B205, 84 (1988). 

[38] C. Lunardini and A. Y. Smirnov, JCAP 0306, 009 (2003). 

[39] C. Lunardini and A. Y. Smirnov, Nucl. Phys. B616, 307 (2001. 

[40] L. Hernquist and V. Springel, Mon. Not. Roy. Astron. Soc. 341, 1253 (2003). 

[41] T. A. Thompson, A. Burrows, and P. A. Pinto, Astrophys. J. 592, 434, (2003). 

[42] R. C. Mitchell et al, Astrophys. J., 574, 293, 2002 

[43] C. Lunardini and A. Y. Smirnov, Astropart. Phys. 21, 703, (2004). 

[44] B. Aharmim et al.,(2005), nucl-ex/0502021. 

[45] T. Araki et al., Phys. Rev. Lett. 94, 081801 (2005). 

[46] Y. Fukuda et al., Phys. Rev. Lett. 81, 1562 (1998). 

[47] Y. Ashie et al., Phys. Rev. Lett. 93, 101801 (2004). 

[48] E. Aliu et al., Phys. Rev. Lett. 94, 081802 (2005). 

[49] R. C. Schirato, G. M. Fuller, (2002), astro-ph/0205390. 

[50] A. Strumia and F. Vissani, Phys. Lett. B 564, 42 (2003). 

[51] B. Jegerlehner, F. Neubig, and G. Raffelt, Phys. Rev. D54, 1194 (1996). 

[52] C. Lunardini and A. Y. Smirnov, Phys. Rev. D63, 073009 (2001. 

[53] H. Minakata and H. Nunokawa, Phys. Lett. B504, 301 (2001). 

[54] M. Kachelriess, R. Tomas and J. W. F. Valle, JHEP 0101, 030 (2001). 

[55] V. Barger, D. Marfatia and B. P. Wood, Phys. Lett. B 532, 19 (2002). 

[56] A. Mirizzi and G. G. Raffelt, |astro-ph/0508612| 

[57] J. M. LoSecco, Phys. Rev. D 39, 1013 (1989). 

[58] T. Yoshida, T. Kajino and D. H. Hartmann, Phys. Rev. Lett. 94, 231101 (2005). 

[59] T. Yoshida, T. Kajino, H. Yokomakura, K. Kimura, A. Takamura and 
D. H. Hartmann, astro-ph/0602195 . 



27 



[60] S. E. Woosley and W. C. Haxton, Nature 334, 45 (1988). 
[61] W. C. Haxton, |a7Xiv:nucl-th/9712049| 
[62] A. M. Hopkins and J. F. Beacom, |astro-ph/0601463 
[63] F. Mannucci et al, Astron. Astrophys. 401, 519 (2003). 

[64] M. Maltoni, T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 6, 122 
(2004) 

[65] J. F. Beacom and M. R. Vagins, Phys. Rev. Lett. 93, 171101 (2004). 

[66] A. M. Dziewonski and S. T. Anderson, Phys. Earth Planet Interiors 25, 297 
(1981). 

[67] M. S. Malek, Ph D Thesis, State University of New York at S tony Brook, Aug. 
2003, available at |http:// www-sk.icrr.u-tokyo.ac.jp/ sk / pub / 

[68] K. Kifonidis, T. Plewa, L. Scheck, H. T. Janka and E. Mueller, 
|astro-ph/0511369 



28 



